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Abstract — We study the clustering properties of advected, non- 
interacting, passive scalar particles in a Burgers fluid witli noise, a prob- 
lem which maps to that of passive sUding particles moving under gravity 
on a surface evolving through the Kardar-Parisi-Zhang equation. Numeri- 
cal simulations show that both the density-density correlation function and 
the single-site mass distribution scale with system size. The scaling func- 
tions diverge at small argument, indicating strong clustering of particles. 
We analytically evaluate the scaling functions for the two-point correlation 
and mass distribution of noninteracting particles in thermal equilibrium in 
a random landscape, and And that the results are remarkably similar to 
those for nonequilibrium advection. 
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I. Introduction 

The advection of a passive scalar in a turbulent fluid is an 
interesting problem which has, in the recent past, attracted con- 
siderable attention (T| f2l f3l. By a passive scalar we mean a 
scalar field (e.g. temperature or density of dye) that is advected 
by the fluid flow but which has no back effect on the flow. In this 
paper we will consider the problem of a passive particle density 
in a compressible fluid described by the Burgers equation in the 
presence of white noise, studied recently in |4|. While an in- 
compressible density field would tend to disperse and ramify 
when carried by an incompressible flow, the reverse can happen 
if the fluid is compressible, and particle clustering (rather than 
dispersal) can result. The characterization of this clustering in 
the one-dimensional case is the principal objective of this paper 

The restriction to flows described by the Burgers equa- 
tion fsl f6l allows the problem to be mapped to another inter- 
esting problem, namely that of passive, sliding particles (called 
sliders) moving under gravity along a growing surface which is 
itself fluctuating in time. Previous studies of the related prob- 
lem of hard-core particles sliding on a fluctuating surface have 
shown the occurrence of a new and interesting state, with large- 
scale clustering of particles [7 | |8|. Now in the problem under 
consideration here, there is no hard-core interaction between 
sliders, so we may expect even stronger clustering properties. 
As discussed below, this expectation is borne out in our numeri- 
cal simulations of the system. There is a clear signature of clus- 
tering in the behaviour of the density-density correlation func- 
tion C(r), which has a scaling form with argument r/L where 
L is the system size. While in the fluctuation-dominated phase 
ordering (FDPO) state discussed in 1 7 1 and 1 8 1, the scaling func- 
tion has a cusp singularity at small argument, we find that it has 
a divergence in the problem under study here. Interestingly, we 
find that this strongly nonequilibrium system has the same scal- 
ing properties as a system of noninteracting particles in thermal 
equilibrium in a random, static landscape. The latter problem is 
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solved analytically, and the results are found to be similar to the 
numerical results for nonequilibrium advection. 

II. Model 

We first discuss the continuum equations describing the pas- 
sive scalar problem with Burgers flow. We then describe a one 
dimensional lattice model which is expected to have the same 
scaling properties. 

The velocity field u of a randomly stirred Burgers fluid is de- 
scribed by the equation 

-^ + \{v-V)v^vV^v + VCh{x,t) (1) 

Here v is the viscosity, while the random stirring is caused 
by Gaussian white noise C,h satisfying {(h{x,t)Ch{x' ,t')) = 
2Dh5'^{x-x')5{t~t'). 

A passive scalar particle moving with the flow follows the 
equation 

dx 

— =av + CS) (2) 
at 

where the white noise C,x{t) represents the randomising effect 
of temperature, and satisfies {Cx{t)Cx{t')) = 2K5{t — t'). The 
parameter a governs the coupling of the particle to the flow. If 
the flow is vortex-free, we may write v = — Vft., where the field 
h satisfies the equation 

^ = iyV'h + ^{Vhf + Ch{x,t) (3) 

This is the Kardar-Parisi-Zhang (KPZ) |9| equation for surface 
evolution with h{x, t) describing the height of a growing surface 
at position x at time t. 

Rather than analysing the coupled Eqs. (2) and (3) directly, 
we will study a lattice model which is expected to have similar 
behaviour at large length and time scales. The model consists 
of a flexible one-dimensional lattice in which particles reside on 
sites, while the links or bonds between successive lattice sites 
are also dynamical variables denoting local slopes. Each link 
takes values +1 (upward slope /) and —1 (downward slope 
\). We evolve the surface by choosing a site at random, and 
if it is on a local hill(^ /\), the local hill is changed to a local 
valley \/). These moves of the links define the well known 
single-step model, whose large distance, long time properties 
are known to be identical to those of the KPZ equation. The 
particles (any number of which can be present at a given site) are 
moved one at a time. We choose a particle at random and move it 
one step downward with probability (1 + K) /2 or upward with 
probability (1 — K)/2. The parameter K, which ranges from 
1 (particles totally following the surface slope) to (particles 
moving independently of the surface) mimic the ratio a/ A in 
Eqs. (2) and (3). In our simulations, we update the surface 
and oarticles at indeoendent sites, reflectina the indeoendence 
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Fig. 1 

The number N(l,L) of particles which lie in the same bin of size l 
as a randomly chosen particle is shown in the inset for various 
system sizes l. the various sizes are 128, 256,512, 1024 as we 

MOVE UP. The MAIN FIGURE SHOWS THE SCALING OF N(L,L) WITH L 



of the noises C,h{^, t) and Cx{t)- This is in contrast to Kardar 
and Drossel's update rule |4| where only particles residing at a 
site affected by the surface evolution are moved. 

III. Numerical Results 

As particles move preferentially downwards, it is evident that 
they tend to move towards valleys, and this valley-seeking ten- 
dency promotes clustering. The question is how one can quan- 
tify its nature and extent. The numerical study of the two- 
point correlation function and single-site particle distribution 
described in this section show that a scaling description affords 
a simple and compact description of clustering properties. 

We consider the steady state of a system of M particles in 
a system with L sites. In our numerical study, we used L in 
the range 128 to 1024, and took M ~ L. The successive up- 
date times of the particle and the surface moves were taken to 
be equal to each other To start with, we monitored the RMS 
displacement xums of a given tagged particle, and verified that 
xrms ~ t^^^ where z = 3/2 is the KPZ value of the dynamic 
exponent. This is expected on the basis that the particle is likely 
to be in the largest fresh valley that has formed within time t in 
the vicinity of the particle, as such valleys have a spatial size 
of order t^/^. Quite direct evidence of particle clustering is 
obtained by choosing a particle at random and monitoring the 
mean number of particles within a distance i of it. The results 
are shown in Fig. 1. Saturation of the curves as £ increases is 
a clear indication of particle clustering. These results are es- 
sentially identical to those of Drossel and Kardar for their 
slightly different model. 

Our results for the two-point (unconnected) density-density 
correlation function G{r, L) — {mimi^r) for different L are 
shown in Fig. 2, where is the number of particles at site i. 
There is strong evidence that the scaling form 

Gir.L) L-'^YirlL) (A) 



is valid with 9 = ^, and that the scaling function Y{y) has a 
power law divergence Y ~ y~'^ as ?/ — > 0. The value of </> 
is close to 1.5, which is the value obtained for particles which 
have reached equilibrium on a random surface, as discussed in 
section IV below. 

The divergence of the scaling function at small argument is 
the result of strong clustering. It is an outcome of the un- 
restricted occupancy of a site by noninteracting particles, and 
quite different from the milder cusp singularity found in the 
FDPO state with hard-core exclusion |8|. As shown in Section 
IV, the form of the scaling function can be found analytically 
for the ostensibly different problem of noninteracting particles 
in thermal equilibrium on a random landscape. As a matter of 
fact, the analytic forms describe quite well the data in Fig. 2 
for the strongly nonequilibrium system under study. This is dis- 
cussed further in Section V. 

Another interesting quantity, which however does not carry 
any spatial information, is the probability P(m, L) that any site 
has occupancy m. As is evident from Fig. 3, P(m, L) follows 
the scaling form 

P(m,L)~-l/(^). (5) 

At first sight, the scaling function seems to follow a power law 
f{y) ^ y"^ with r close to 1, but guided by the analytic theory 
to be discussed in the next section, we find that the ?/ — > be- 
haviour is well described by the form f{y) ~ \n{b/y) with 
6 ~ 5. On integrating over m, we find 1 — P(0, L) ^ ■i(lni)^, 
implying that the number of occupied sites increases with size as 
(InL)^. This indicates that the particles are distributed within a 
relatively compact region in space. Finally, we note that Eq.(5) 
leads to (m^) = G{0, L) ^ L. Direct simulation results for 
G(0, L) verify the linear dependence on L, lending further sup- 
port to the scaling form (Eq.(5)) 

IV. Analytical Results (Equilibrium) 

Since the particles are noninteracting in our model, the num- 
ber density p{r,t) — m{r,t)/L of particles is identical to the 
probability that a single particle, moving in the fluctuating land- 
scape, will be at a position < r < i at time t. Thus it suf- 
fices to consider one single particle moving in the fluctuating 
landscape. We consider the limiting case when changes in the 
surface occur over a time scale much larger than the time scale 
of fluctuations of the particle. In this limit, for any given surface 
configuration h{r), the Langevin equation (Eq.(2)) describes the 
evolution of a particle in a potential h{r). Before the surface 
configuration h{r) can change, the particle will relax to its its 
equilibrium Gibbs state, where the probability of finding the 
particle at position r is given by p{r) = cxp[— /3/i(r)]/Z with 
the partition function Z = exp[~(3h{r')]dr' , where /3 = 
l/ksT. Clearly p{r) is a random variable which changes from 
one configuration of h{r) to another. We assume that the surface 
itself has reached the stationary state of Eq. (3). It is well known 
that for the 1-d KPZ equation, the stationary state is described by 
the following measure, Prob[{ /i(r}] cx exp [— ^ / h'^{r')dr''\. 
Thus, any stationary configuration can be thought of as the 
trace of a random walker in space evolving via the equation, 
dh(r) Idr = £ (r) where the white noise £ (r) has zero mean and 
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is delta correlated, {^{r)£,{r')) — 6{r — r'). Thus our problem 
then reduces precisely to the celebrated Sinai model 1 1 1 1 where 
a particle moves in a random potential which itself is a random 
walk in space. There is a slight difference however. The periodic 
boundary condition on the surface, h{r) — h{r + L), indicates 
that we are considering a random potential which is pinned at its 
two ends, i.e. a Brownian bridge, rather than a free walk. 

Given this stationary measure of the surface configurations, 
we then want to evaluate the correlation function {p{ro)p{r + 
To)) = G{r, L)/ L"^ (evidently independent of tq due to trans- 
lational invariance) where the angular brackets denote the aver- 
age over the surface configurations sampled from the stationary 
measure mentioned above. Thus we have 



L-^G{r,L) = 



.-Pmro)+h(ro+r)] 
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(6) 



Fortunately, the object on the right hand side of Eq.(|6} was eval- 
uated exactly by Comtet and Texier 1 10| in the completely dif- 
ferent context of one dimensional disordered supersymmetric 
quantum mechanics, where the right hand side of Eq. is sim- 
ply the correlation function in the ground state wave function. 
They actualy evaluated the n-point correlation function. Adapt- 
ing their results to our case with n = 2, we find the following 
expression 



r/L 
Fig. 2 

Scaling of the density-density correlation function G(r,L) 

WITH system size L. THE DASHED CURVE IS A PLOT OF THE SCALING 
FUNCTION Y (EQ.(8)) SHIFTED DOWNWARDS FOR CLARITY. 



G(y) = —K,{y). 
V 



(10) 



Note that the expression for G(r, L) has the expected symmetry, 
G(r, L) = G{L — r, L). It is easy to evaluate G(0, L) from Eq. 
(|8j for large L. We find, G(0, L) w l3^L/12 for large L. On the 
other hand, in the scaling limit, r ^ oo, L ^ oo but keeping the 
ratio y — r/L fixed, one finds from Eq. (7) that for < a; < 1, 
G(r, L) ~ L~^/'^Y(r/L) where the scaling function is given by 



(8) 



Note that the point r = is not part of the scaling function. 
In fact, the power law divergence of the scaling function as 
y ^ 0, Y{y) ~ with (j) = 3/2, is necessary in order that 
G(r ^ 0, L) - L for large L. 

This formalism can also be used to calculate the equilibrium 
probability density P{p, L). Details of the calculation will be 
given elsewhere. Our results indicate that P(/0, L) which can be 
written as the sum of two parts : 



(9) 

The first part refers to vacant stretches, and to the fact that 
the number of occupied sites occupies a vanishing fraction 
^ (lni)^/i of the system. The scaling function G(j/) in the 
second oart is eiven bv 



^-^yhere Ko{y) is the modified Bessel function which has the 
asymptotic behaviour [— ln(?;/2) — 0.5772...] as y ^ 0. The 
theta function incorporates a lower cutoff on the validity of the 
scaling form. 

The equilibrium-based results obtained in this section appear 
to describe rather well the behaviour of the nonequilibrium sys- 
tem studied numerically in section III. For instance Eq.(8) de- 
scribes the behaviour (Eq.(4)) of the two-point correlation func- 
tion (Fig. 2), with /3 ~ 4 giving good agreement. The data for 
the mass distribution (Fig. 3) is quite well described by Eq.(lO), 
but with a different value of /3 (~ 2.3). An intermediate value 
of /3(~ 2.7) gives tolerable, but not very satisfactory, agreement 
with both sets of data (Figs. 2 and 3). 

V. Discussion 

As we have seen in the previous two sections, a number of 
the properties of the nonequilibrium system of particles sliding 
in a fluctuating landscape are strikingly similar to those of an 
equilibrium system of particles settling in a quenched random 
landscape. This raises two questions: First, what is the source 
of this similarity in behaviour? And second, how robust is it? 
We address these below. 

In thermal equilibrium, in a given configuration of the land- 
scape, the Boltzmann factor guarantees that most particles are 
to be found within a height T of the global minimum, where 
T is the temperature. In the nonequilibrium situation, by con- 
trast, examination of several typical configurations shows that 
particles cluster near valley bottoms, but these valleys are of- 
ten not the globally lowest one. Why then are properties simi- 
lar when averaged over configurations? A oartial answer lies in 
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The PROBABILITY P(M,L) OF FINDING A SITE CONTAINING M PARTICLES 
IS SHOWN IN THE INSET FOR VARIOUS SYSTEM SIZES L. THE MAIN FIGURE 
SHOWS THE SCALING OF P(M,L) WITH L; THE DASHED CURVE SHOWS THE 
SCALING FUNCTION G (EQ.(IO)) SHIFTED DOWNWARDS FOR CLARITY. 



the fact that terrain around any reasonably deep valley is statisti- 
cally similar. This feature of the 'random walk' landscape under 
consideration can plausibly lead to particles being distributed in 
similar fashions in the two cases. 

Evidently, this happens despite the fact that the character of 
the noise which leads to particle re-settling is quite different in 
the two cases: in the equilibrium problem, particles are acted 
upon by thermal noise, while in the nonequilibrium case, it is 
surface fluctuations which drive re-settling. Thus, despite the 
similarities in the landscape, it is not clear to what extent one 
should expect an equivalence of results, and indeed, as the com- 
parison in Section IV shows, the equivalence for different prop- 
erties involves equilibrium systems at different temperatures. It 
would also be interesting to see how robust these results are with 
respect to variations of the parameter oj, which is the ratio of the 
times Ts between successive updates of the surface, and Tp be- 
tween successive updates of the particles. In the studies reported 
in this paper, lu was held constant at 1. Preliminary investiga- 
tions suggest that variations of oj can induce interesting features 
specific to the nonequiUbrium problem. 
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